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ABSTRACT 

In this paper we model the gravitational wave emission of a freely precessing neu- 
tron star. The aim is to estimate likely source strengths, as a guide for gravitational 
wave astronomers searching for such signals. We model the star as a partly elastic, 
partly fluid body with quadrupolar deformations of its moment of inertia tensor. The 
angular amplitude of the free precession is limited by the finite breaking strain of the 
star's crust. The effect of internal dissipation on the star is important, with the preces- 
sion angle being rapidly damped in the case of a star with an oblate deformation. We 
then go on to study detailed scenarios where free precession is created and/or main- 
tained by some astrophysical mechanism. We consider the effects of accretion torques, 
electromagnetic torques, glitches and stellar encounters. We find that the mechanisms 
considered are either too weak to lead to a signal detectable by an Advanced LIGO 
interferometer, or occur too infrequently to give a reasonable event rate. We therefore 
conclude that, using our stellar model at least, free precession is not a good candidate 
for detection by the forthcoming laser interferometers. 

Key words: accretion - radiation mechanisms: non-thermal - relativity - 
starstmagnetic fields - stars: neutron - stars: rotation 



1 INTRODUCTION 



Freely precessing neutron stars have long been recognis ed as 



a potential sou r ce of detectable gra vitational waves (Zim- 
[mcrmann 1978 ; Alpar fc Pines 1985). Des p ite their regula r 
inclusion in review articles (Thorne 1987; Flanagan 199^) 
they have received little in the way of detailed modelling. In 
this paper we will combine many of the physical processes 
relevant to the problem of free precession and assess their 
relative importance. The motivation b ehind this w ork is to 
aid gravitational wave data analysis (Schutz 1991). Given 
the huge computational requirements of this analysis, any 
additional information supplied by theoretical modelling of 
a source greatly increases the chances of detection. With the 
TAMA detector already operational, the GEO600 and LIGO 
detectors due to go on-line within a year, and the VIRGO 
detector following soon after, this issue is particularly press- 
ing. 

The study of gravitational wave generation from freely 
precessing neutron stars can be divided into two parts. The 
first problem is the formulation of a free precession model 
consistent with our current understanding of neutron star 
structure. Neutron stars are not rigid bodies — they consist 
of a thin elastic shell containing a superfluid core. A model 
taking these features into account was described in detail 
in Jones & Andersson (2001), where the effect of free pre- 



cession on electromagnetic pulsar signals was described and 
compared with pulsar observations. (A handful of potential 
free precession candidates were identified, but, as we will 
see, they all rotate too slowly to be of gravitational wave 
interest). A brief summary of this model is set out in this 
paper. A key feature is the decay of free precession due to 
dissipative processes internal to the star. 

The second part of the problem is to look at particular 
scenarios in which free precession is created and/or main- 
tained by an astrophysical mechanism. The torques due to 
accretion disks, neutron star magnetic dipole moments, and 
other gravitating bodies will be considered, as well as per- 
turbations associated with glitches. We wish to investigate 
whether these mechanisms are capable of balancing the in- 
ternal dissipation to give steady long-lived precessional mo- 
tions. Having done so, we will then be in a position to esti- 
mate possible gravitational wave amplitudes in these various 
scenarios. 

The structure of this paper is as follows. In section ^ we 
briefly describe our free precession model, and parameterise 
the timescale in which internal dissipation damps the wobble 
motion. In section ^ we show how the finite breaking strain 
of the crust can be used to place an upper bound on the 
gravitational wave field of a precessing star, regardless of 
its environment. We also suggest a detection strategy for 
gravitational wave data analysts. General formulae relating 
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the gravitational wave amplitude to the strength and nature 
of a pumping torque are given in section ^ The possible 
gravitational wave field strengths in various scenarios are 
then estimated in sections BMsI Our conclusions are given in 
section ^, together with some suggestions for further work. 



2 A MODEL OF NEUTRON STAR FREE 
PRECESSION 

2.1 Dynamics of free precession 

In this section we will briefiy summarise our model of neu- 
tron star free precession. For more detail see Jones & Ander- 
sson (2001). We will begin by describing the free precession 
of a rigid body, as the motion in the more realistic elastic 
shell/fiuid core case can be thought of as a modification of 
this. 

The moment of inertia tensor of any axisymmetric rigid 
body can be written as 



I = /o<5 + A/d(ndnd-5/3), 



(1) 



where 6 is the Kronecker delta, and the unit vector nj points 
along the body's symmetry axis. Then the principal mo- 
ments are h = h = lo — AJd/3, I3 = lo + 2A/d/3, so that 
I3 — Ii = A/d. When A7d > the body is said to be oblate, 
and when negative the body is prolate. As we will describe 
below, the oblate case is the more physically plausible. 

The angular momentum is related to the angular veloc- 
ity according to 



J = (Jo - A/d/3)fl - A/dJ^and, 



(2) 



where the 3-axis lies along nd. This shows that the three 
vectors J, CI and nd are always coplanar. Following Pines 
& Shaham (1972a,b) we will call the plane so defined the 
reference plane (see figure ^). Given that the angular mo- 
mentum is fixed, this plane must revolve around J. The free 
precession is conveniently parameterised by the angle 9 be- 
tween nd and J. We will refer to this as the wobble angle. 
For a nearly spherical body the angle 9 between fl and J is 
much smaller than the angle between J and nd, according 
to 

a ~ — — sm o cos 9. 

Ii 



(3) 



We will denote by nj the unit vector along J. Decomposing 
the angular velocity according to 



fl = (pnj + 7/>nd 
then gives 



J = /i 



A/d 



(4) 

(5) 
(6) 



The symmetry axis nd performs a rotation about J in a 
cone of half-angle 9 at the angular frequency tp. We will 
refer to this as the inertial precession frequency. There is 
a superimposed rotation about the symmetry axis nd at 
the angular velocity ^ . This is usually referred to as the 
body frame precessional frequency, with the corresponding 
periodicity known as the free precession period: 

P^. - ^. (7) 




Figure 1. This figure shows the reference plane for a freely pre- 
cessing body, which contains the deformation axis nd , the angular 
velocity vector f2 and the fixed angular momentum J. The vectors 
nd and f2 rotate around J at the inertial precession frequency rp. 
We refer to B as the wobble angle. The vector Txj is the part of 
an applied torque that causes a secular evolution in the wobble 
angle. 



For a nearly spherical body equation (g) shows that ip 4>- 
Note that the angles {9, (p, ip) are simply the usual Euler 
angles which describe the orientation of the rigid body (see 
e.g. Landau & Lifshitz 1976, figure 47). 

Turning now to the more realistic case of an elastic crust 
with a liquid core, the moment of inerti a of the crustal shell 
can be written as (Alpar & Pines 1985): 

I^IoS + A/n(nnnn - S/3) + AJd(ndnd - S/3). (8) 



The first term on the right hand side is the moment of inertia 
of the non-rotating undeformed spherical shell. The second 
term is the change due to centrifugal forces, and has no, the 
urnt vector along CI, as its symmetry axis. The third term is 
the change due to some other source of deformation, (such 
as strains in the crustal lattice), and has the unit vector nd, 
fixed in the crust, as its symmetry axis. Without this third 
part the above equation would simply represent a fiuid ball, 
and free precession would not be possible. 

The free precession of such an elastic shell containing 
a liquid core is then very similar to the rigid result, with 
the geometry described in figure 1 still applying, providing 
the A/d in equations (^ and (^ is now set equal to the 
deformation in the inertia of the whole star, i.e. crust and 
core, while /o is still equal to the crustal moment of inertia 
only. Explicitly: 

A/d , 



T ' 

-» crust 

A/d 



(9) 



(10) 



/crust 

Note that centrifugal deformation A/n is not of importance 
when considering the free precession geometry. 

One further component can be added to our model: A 
pinned superfiuid coexisting with the inner crust. As was 
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first described by Slialiam (1977), tlie effect of sucli a compo- 
nent is to increase tlie body frame free precession frequency 
tp. In tlie case where tfie rotation rates of tfie superfluid and 
crust are tlie same, equations and (^o|) still apply, with 
Aid now containing a part equal to the moment of inertia 
of the pinned superfluid. However, as described in Jones & 
Andersson (2001), the handful of free precession candidates 
identified in the pulsar population don't seem to have such a 
component. Therefore, on those few occasions in this paper 
where we assume a particular source of deformation, we will 
assume that Aid is caused entirely by Coulomb forces in 
the crustal lattice. (If real stars do have a significant pinned 
superfluid component, then the increased body frame pre- 
cession frequency will lead to an even faster dissipation of 
the free precession energy, and tend to re-enforce our final 
conclusion) . 

The effect of a torque on the free precession can be 
easily summarised. Given that the torque-free motion is de- 
termined completely by the two numbers {4>,9), we need 
only describe the effect of the torque on these. If the torque 
causes the magnitude of the angular momentum to change 
at a rate J then 

'<i> = r- (11) 
Ji 

The evolution in the wobble angle is determined by the com- 
ponent of the torque projected into the reference plane which 
lies perpendicular to J: 

e = (12) 

J crust 

as illustrated in figure m. 



energy of the star. The constant B also depends on the equa- 
tion of state, and will be of order of the total electrostatic 
binding energy of the ionic crustal lattice. We will define 
b = B/{A-\- B) as the rigidity parameter. It is equal to zero 
for a fluid star [B — 0) and unity for a perfectly rigid one 
{B/A — » oo). Realistic neutron star equations of state imply 
that b takes a value of: 

1.6 X 10''^ rI/mI^. (15) 

(See Jones (2000) for a simple derivation, and Ushomirsky, 
Cutler & Bildsten (2000) for a detailed numerical treat- 
ment). The smallness of this number reflects the fact that 
gravitational forces dominate crustal Coulomb forces in de- 
termining the equilibrium shape of the star. In this sense, 
neutron star crusts are very far from perfectly rigid. Note 
that b ~ 10~^ for a canonical 1.4Mq, 10 km neutron star. 
More rigid stars can exist only if less massive (and there- 
fore larger radii) neutron stars occur in nature, or if current 
equations of state seriously underestimate the crust thick- 
ness. We will therefore present results for rigidity parameters 
over the interval 10""^ 10~^, but will bear in mind that 
values at the smaller end of the interval are more plausible. 

For a very young star the zero-strain oblateness eo will 
be determined simply by the star's shape at the moment 
when its crust first solidified, i.e. eo ~ en(/ ~ /solid), 
where /soud is the spin frequency at the moment of so- 
lidification. For older stars, eo will have changed due to 
plastic deformation in the crust, either because of a grad- 
ual creep or a more violent shape change, possibly con- 
nected with a glitch. Certainly, it seems likely that real stars 
will be oblate (eo > Aid > 0) rather than prolate 
(eo < 0^ A7d < 0). 



2.2 Sources of deformation 

The centrifugal deformation described above can be con- 
veniently parameterised by the dimensionless quantity en 
which we will define by 3en/2 = A/n//star (this follows the 
notation of Alpar & Pines 1985). This will be of the order 
of the rotational kinetic energy of the star divided by its 
gravitational binding energy: 



eo 



GM 



2.1 X 10" 



(loik)^^^/^^ 



(13) 



where Q = 2nf is the angular frequency, i?6 the neutron 
star radius in units of lO^cm, and Mi. 4 the mass in units of 
I.AMq. 

The deformation Aid we can similarly parametrise in 
dimensionless form using the relation 3ed/2 = AJd/^star. 
This deformation is due to some physical process other than 
rotation, which need not be specified in many of our gravi- 
tational wave estimates. However, in practice the most sig- 
nificant source of deformation in a neutron star is likely to 
be strains in its solid crust. As described in Baym & Pines 
(1971) and Pines & Shaham (1972a,b), if the crust has a 
zero strain oblateness eo, the actual deformation produced 
is of order 



B 



ed 



A + B 



eo- 



(14) 



where the constant A depends on the stellar equation of 
state, and will be of the order of the gravitational binding 



2.3 Wobble damping 

A real neutron star, once set into free precession, will not 
precess forever — energy will be dissipated within the star, 
converting the kinetic energy of the wobble into thermal en- 
ergy. Also, gravitational wave energy and angular momen- 
tum will be radiated to infinity, which must be subtracted 
off the star's motion. 

The problem of gravitational radiation reaction was ex- 
amined in detail by Cutler & Jones (2000), who showed that 
the main result was to cause the wobble angle to decay ex- 
ponentially on a timescale: 

10 1 n9 / ^crust 

Tg = 1.8 X 10 yr ' 



V 10'*''gcm2 



lO^^gcm^ 

aL 



/ lOOHz 



(16) 



regardless of whether the deformation is oblate or prolate. 

This is almost certainly much longer than the timescales 
connected with internal dissipation. In particular, models 
of neutron star interiors, motivated in part by the need to 
model glitches, predict a frictional type coupling between 
the crust and core. The result of this coupling is a torque T 
exerted on the crust: 

T = A'(nfluid - fisoiid), (17) 
where K is a, positive constant. Such a torque would tend 
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to restore corotation between crust and core in a glitching 
neutron star. However, as modelled by Bondi & Gold (1955), 
the torque would also tend to damp the wobble angle of 
a precessing body. The timescale for this damping can be 
parameterised by n, the number of free precession periods 
Pfp in which one e-fold occurs: 

^°-nP. (18) 



Td = 



where P is the spin frequency (approximately 2-K/ip) of the 
body. The parameter n has been estimated by Alpar & Sauls 
(1988) who examined the scattering of electrons off the su- 
perfluid vortices. This interaction is sometimes referred to 
as 'mutual friction'. Alpar & Sauls estimated n to lie in the 
interval 400 10'*, giving a wobble damping timescale of: 



Td 



3.2 yr 



(loO 



100 Hz 



/ 



10*4 gcm2 



10^* g cm^ 



(19) 



Comparing equations ( p^ ) and (^^ it is clear that internal 
dissipation is likely to damp the free precession much more 
rapidly that gravitational radiation reaction. 

Interestingly, in the case of prolate deformations, the 
internal dissipation acts so as to i ncr ease the wobble angle. 



However, as described in section 2.2, prolate deformations 



are probably not likely to be found in real stars, and so for 
the remainder of this paper we will assume that deforma- 
tions are oblate, corresponding to a damping of the wobble 
motion. 



3 THE GRAVITATIONAL WAVE FIELD 
3.1 General form of the field 

The gravitational wave field for a rigid precessing body was 
first calculated by Zimmerman & Szedenits Jr (1979), us- 
ing the mass quadrupole formalism (see e.g. Misner, Thorne 
and Wheeler 1973). As described in the previous section, 
even though the elasticity and fluid core are important in 
determining the free precession period, the geometry of free 
precession in the realistic case is very similar to that in 
the rigid case — the deformation bulge moves in a cone of 
half-angle S at a rate 0. (The superimposed rotation at 
about the axis rid does not change the mass quadrupole of 
the body, and so does not appear in the gravitational wave 
physics). It follows that, to the accuracy of our model at 
least, the wave field calculated by Zimmerman & Szedenits 
Jr (1979) applies in the realistic case. Explicitly, the gravi- 
tational waves are emitted at frequencies cj) and 20 and with 
the two polarisations (denoted by + and x): 



h+{4>) 



2(jf 



sin i cos iA/d sin 6 cos 6 cos(<^t) 



hx {4>) = sin iA/d sin cos 9 sm((l>t) 



h+{2(j>) = — (H-cos^ i)A/dSin^6lcos(2(^f) 
r 

hx(2<j}) = ^2cosiA/dSin^6lsin(2<^t) 
r 



(20) 
(21) 
(22) 
(23) 



for a source at distance r and angular momentum at inclina- 
tion angle i to the line of sight. Note that the A7d factor is 
the part of the quadrupole moment tensor due to the defor- 
mation process, not the centrifugal piece. The modifications 
to this result due to the non-rigidity are small, and lie be- 
yond the accuracy of our free precession model. 

3.2 Limit on gravitational wave amplitudes due to 
finite crust strength 

As described in Jones & Andersson (2001), a real neutron 
star crust will have a finite breaking strain iibrcak- This can 
be used to place an upper bound on the wobble angle of 
a precessing star. We know that (for small wobble angles 
at least) when an initially axially symmetric body is set 
into free precession, a deformation A/d remains along the 
axis Ud fixed in the star, while a deformation AIq points 
along the angular velocity vector. From the point of view 
of an observer attached to the crust, a deformation of size 
AIq describes a cone of half-angle 6 + 9 ~ 9 about Ud- 
This change in shape is all we need to know to estimate 
the strain: The change in position of any given particle is of 
order Ren9, while the corresponding strain is of order en9. 
This precession-induced strain is not constant, but varies 
with magnitude enO over one (body frame) free precession 
period. As there exists a maximum strain Mbrcak that the 
solid can withstand prior to fracture, the wobble angle will 
be limited to a value of Ubroak/en so that: 



; 0.45 



100 Hz 



'^brcak 

10-^ 



radians. 



(24) 



Qualitatively, we can say that faster spinning neutron stars 
have larger bulges to re-orientate and therefore can sus- 
tain smaller wobble angles prior to fracture. For sufficiently 
slowly spinning stars the above equation breaks down, yield- 
ing angles in excess of 7r/2. The wobble angles of such slowly 
spinning stars are not limited by crustal strain. 

The value of Ubroak is highly uncertain. By extrapolat- 
ing the breaking strains of terrestrial materials, Ruderman 
(1992) suggests the value relevant for neutron star crusts 
may lie in the range 10~^ to 10"'*. 

This limit on the wobble angle 9 can be used to place 
an upper bound on the gravitational wave amplitude of a 
freely precessing star with a given deformation A/d. To do 
so we will characterise the field strength by 



h : 



G 
7^' 



-A/d 



(25) 



This is an order-of-magnitude approximation to the set of 
equations (|2^)-(^). The trigonometric factor describing the 
variation of wave amplitude with inclination angle i has been 
neglected entirely, while the trigonometric factor in the pre- 
cession angle 9 has been replaced by its small angle limit. 
Note that in this small-angle limit the radiation is emitted 
mainly at the frequency (p — the radiation at frequency 20 is 
a factor 9 smaller. 

The purpose of this approximation is to enable us to 
gain a insight as to how gravitational wave amplitudes de- 
pend upon source parameters such as the breaking strain, 
frequency and deformation size. This analysis is a useful 
preparation for sections where detailed astrophysical sit- 
uations are considered. 
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We will specialise to the case where the deformation 
A7d is due to Coulomb forces in the crustal lattice. We set 
A7d = 3/stared/2 with Ed given by equation (p^. For defi- 
niteness, we will assume the crust is 'relaxed', i.e. the zero- 
strain oblateness eo is equal to the centrifugal oblateness en . 
(These two quantities can differ by the small crustal break- 
ing strain anyway — see Jones & Andersson 2001). We then 
find 

»-'"(io^)'(t^)(^)» '»> 

for a spin frequency /. In ( p6|) has been left a free pa- 
rameter. However, to obtain an upper bound on h we can 
set 9 to its maximum value as obtained from crust cracking 
considerations, i.e. put 6 = Smax. We then obtain 

M^ = ^.a.)=3.7xl0-(^)'(^) (27) 
for / < fe and 

M»^»,„„)^i.8x,o-»(^)'(jL,)(iipf),a, 

for f > fe, where fg is the frequency at which ^max ~ 1: 
/«-69Hz(^)'^' (29) 

For spin frequencies close to and below this our small-angle 
approximation breaks down, so we have put 6 = 1 for f < fe 
so as to still obtain results correct to within an order-of- 
magnitude. 

In figures ^ and ^ we have plotted the maximum ampli- 
tude h for a variety of neutron star parameters. The noise- 
curves have been taken from Owen & Sathyaprakash (1999). 
Note that a knee appears in many of the signal curves. This 
knee corresponds to f = fg. Above this frequency the wob- 
ble angle is limited according to (^^. We have assumed the 
matched filtering can accumulate signal only for a time of 
one year. The amplitudes are shown only for frequencies 
less than a kilohertz, as stars rotating more rapidly than 
this have not yet been observed. Figure ^ plots h for the 
parameters Ubreak = 10~^, r = 1 kpc and for three dif- 
ferent values of 6: 6 = 10"^ 6 = 10"*, b = 10'^. RecaU 
that b = 10~^ is the value expected for a canonical IAMq 
R = 10 km neutron star, while b = 10~^ would correspond 
to a lighter larger star with a thicker crust. In figure M plots 
h for 6 = 10~^ and breaking strains of Mbreak = 10 , 10~^ 
and 10"''. 

In making the above assumptions regarding matched 
filtering we have assumed that the wobble angle and spin 
frequency of the star remain constant for the observation 
period of one year. In practice the wobble an gle will decay 
significantly over this interval (see section |2.3| ) unless a suf- 
ficiently strong pumping mechanism is active. However, the 
purpose of these figures isn't to model in detail any particu- 
lar scenario — they are included simply so that we might gain 
insight into how wave amplitudes depend upon the rigid- 
ity and breaking strain parameters, and identify parameters 
values capable of leading to detectable signals. These figures 
say nothing about the distance to the nearest source, or, in 
the case of burst-like sources, the event rate. Issues such as 
these can only be addressed in the case of particular pump- 
ing mechanisms. This will be carried out in sections H-fel 
where the complicating effect of internal dissipation will be 




12 3 4 

log(f) 



Figure 2. The maximum gravitational wave amplitude for 
Coulomb deformations. The star is at a distance of 1 kpc with 
"break = lO'^ and b = 10-3, 10"*, 10-5. The matched filter has 
been assumed to accumulate signal for an interval of one year. 
The noisecurves have been taken from Owen & Sathyaprakash 
(1999). 




log(f) 

Figure 3. The maximum gravitational wave amplitude for 
Coulomb deformations. The star is at a distance of 1 kpc with 
b = 10"5 and Mbroiik = 10"^, 10" 10" The noise curves are as 
indicated in figure ti. 

included also. In short, figures |^ and ^ represent absolute 
upper bounds on the wave amplitude for a given spin fre- 
quency, crustal breaking strain and deformation, but give 
no information about likely wave amplitudes in nature. 

With these qualifications in mind, we will note that 
these figures show that for neutron stars with b ^ 10"'', 
the 9 = Smax gravitational wave signal is detectable by first- 
generation interferometers for frequencies in excess of 100 
Hz. Slower spinning sources produce steady gravitational 
wave signals which are intrinsically too weak to be detected. 
Stars with 6 ~ 10"* are potentially detectable by first gen- 
eration interferometers only for frequencies close a 1 kHz. 
More realistically, stars with b ~ 10"^ are potentially de- 
tectable by an Advanced LIGO interferometer at frequencies 
over 100 Hz. 
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3.3 Detection strategy 

The problem of determining source parameters from elec- 
tromagnetic observations of a freely precessing neutron star 
was examined by Jones & Andersson (2001). The extrac- 
tion of these parameters in the case of a gravitational wave 
detection was examined by Zimmermann & Szedenits Jr. 
(1979), who showed that measurement of all four compo- 
nents of h given in equations (|20[)-(|2^) allows extraction of 
the inclination angle i, the wobble angle 9 and the quantity 
A7d/r. 

As was also noted by Zimmermann & Szedenits Jr. 
(1979), if this is combined with an estimate of r a value for 
A7d is obtained. This is of considerable use: If the deforma- 
tion is of Coulomb origin the equation Aid ~ 3/star6eo/2 can 
be combined with equation-of-state calculations of the rigid- 
ity parameter b to give the reference oblateness eo. A value 
of eo very different from the fluid oblateness might then be a 
sign that the deformation Aid is not due to crustal strains. 

Note however that if the wobble angle is very small the 
above two equations can be solved only for i and the product 
AldO/r. This is likely to be the case for the fastest spinning 
stars whose wobble angles are limited by fracture according 
to equation (p^. 

If the source is observed as a pulsar it may be possible to 
obtain additional information. Suppose the electromagnetic 
pulse originates from a dipole m inclined at an angle x to 
the deformation axis. Then the motion of m is due to the 
combined rotation of rid about J at and the motion of m 
about Ud at tp. As noted by Zimmermann & Szedenits Jr. 
(1979) the electromagnetic pulse frequency is then given by 
/cm — (j> + ip tor 9 < X, while for > x we have /cm = 4>- 

In the case of Coulomb deformations of rapidly spinning 
stars we would expect the 9 < X regime to apply. Then the 
difference between electromagnetic and gravitational wave 
frequencies would give ip, and to leading order equation (jlo|) 
gives 



A7d = /crust V'//ci 



(30) 



If in addition an estimate of the pulsar's distance were avail- 
able the product {AId/Iciust)9/r could be decomposed into 
its component parts. 

The above remarks concerning frequency splitting sug- 
gest the following strategy for gravitational wave observers. 
Suppose we are searching for a precessional gravitational 
wave signal from a star of known (average) electromagnetic 
pulse frequency /. Then one matched filter should be used 
at that frequency to cover the 9 > x case where there is 
no splitting between the gravitational and electromagnetic 
signals. However, if splitting does occur, it will be necessary 
to search at a rotational frequency A/ = /A/d //crust above 
this. For definiteness we will assume a Coulomb deformation 
to give 



A/ 



(nil 



10 



(31) 



where equation ( p^ ) has been used. 

We now wish to identify a frequency band in which 
to search for the gravitational signal. This is easily done. 
Consider some particular frequency, say 100 Hz. Then from 
figure ^ we see that even for an Advanced LIGO, b values 
less than 10~* give signals too weak to be detected. From the 



above formula this value of b corresponds to A/ = 3 x 10~^ 
Hz. Therefore there is no point in searching for signals at 
frequencies less than (100 -I- 3 x 10~^) Hz (apart from the 
signal at / mentioned above). On the other hand, on physical 
grounds the value b = is surely an upper bound on 

the rigidity parameter. This corresponds to A/ = 3 x 10~^ 
Hz. Therefore there is no point in searching for signals at 
frequencies greater than (100 -I- 3 x 10~^) Hz. In this way 
we identify a frequency band in which to fruitfully search 
for the gravitational signal. A year-long integration would 
require matched filters with a 3 x 10~* Hz spacing. In this 
example this would require approximately 10"" templates. 

The above argument attempts to identify the gravita- 
tional waves at the inertial precession frequency (j). There 
will be radiation at 20 also. This can be searched for using 
templates at twice the frequencies described above. 



4 CLASSIFICATION OF PUMPING 
MECHANISMS 

It is useful to categorise pumping mechanisms according to 
the sort of evolution they produce in the wobble angle. As 
discussed above, the effect of a torque depends upon its pro- 
jection into the reference plane, which rotates around the 
star's angular momentum vector at a rate 0. This projection 
will depend upon the details of the mechanism producing the 
torque, possibly leading to a complicated evolution in the 
wobble angle. Nevertheless, we will find it useful to define 
the four categories described below. We will define Tpump 
as a timescale characterising the strength of the pumping 
torque. Its exact meaning will depend upon the nature of 
the mechanism. 

It is also useful to derive the following expressions for 
the wave amplitude written in terms of a timescale rather 
than a frequency. Equation ( |25| ) can be combined with equa- 
tion (|l|) to give the wave amplitude in terms of the gravi- 
tational alignment timescale: 



h = 



-^crust 

2c3 



Tg.e 



1/2 



(32) 



This will be extremely useful when estimating wave ampli- 
tudes for particular pumping mechanisms. Combining equa- 
tion (js^) with equation (^J) then gives an estimate of the 
maximum possible wave amplitude at a given frequency and 
gravitational alignment timescale Tg^g: 



= 1.4 X 10 



-24 



lO'^'gcm^ 



1/2 



/ 10^yr\^''^ /lkpc\ 
l^j (^) 



(33) 



for f < fe, and 



/imax = 6.7 X 10 



25 / -t crust 



10«gcm2 



1/2 



/ ^tbroak \ 



100 Hz \ / lO'' yr \ ' fl kpc 



1/2 



(34) 



for f>fe 
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In a similar way equations ( p5| ) and can be com- 
bined to give the wave amplitude in terms of the internal 
dissipation timescale ra: 



2ttG n(l) e 

— Jcrust • 

C* Td r 

Setting 9 equal to its maximum value then gives 



(35) 



/imax = 3.3 X 10 



-30 / 



(lOO Hz) 



/ /crust \ / 10^yr\ /Ikpc 



(^) (^) (»> 



yiO«gcm2 
for f < fe, and 

f, 1 R ^ in-30 f ^ \ / -"break \ / 100 Hz 

/..ax = 1.6X10 

-^crust 



10*4 g cnr 



for f > fe- Again, these formulae will be of use in estimating 
wave amplitudes in later sections. 



4.1 Oscillatory pumping 

We will first consider the case where the torque does not 
remain fixed relative to the reference plane, but instead ro- 
tates around it at some frequency. In sections ^ and ^ we 
will identify such torques that rotate at the body frame pre- 
cessional frequency with respect to the reference plane, i.e. 
with a frequency ip. Specialising to this case we see that for 
half a precession period the wobble angle will increase, but 
for the next half-cycle it will decrease. We will refer to this 
as oscillatory pumping. It will produce a wobble angle that 
varies as 



6 = doscil + costpt) 



(38) 



where Ooac denotes the average value of 8 attained. To order- 
of-magnitude accuracy this will be given by 

r, /^fp /crust 1 



(39) 



where equation ( |lO| ) has been used. The wave field is then 
estimated by substituting ( ^jj ) into ( pj] ) to give 



h ' 



2ttG Jc, 



(40) 



Note that this is independent of the size of the deformation 
Aid- However, in the case of sufficiently small Aid the pre- 
cession angle as g iven by ( ^9| ) will exceed Smax- In such a 
case equation (WOl) represents an upper bound. 



4.2 Exponential pumping 

In subsequent sections we will find examples of torques 
which tend to increase the wobble angle at a rate propor- 
tional to the angle itself, i.e. produce exponential pumping. 
Then 



1 



1 

Td 



(41) 



In general, these timescales are all functions of frequency. In 
the case where the frequency and therefore the timescales 
are not constant the solution will be complex. However, in 
a constant frequency system such as an accreting star at 
spin equilibrium, the above equation leads to exponential 
solutions. We will consider such a system. 

Clearly, in order to calculate the wave amplitude due to 
such a mechanism we need to compare the three timescales 
that appear in (pl|). In order to simplify this problem we will 
divide it into two smaller ones. First we will include only the 
pumping and gravitational radiation reaction. This will lead 
to an upper bound on h for stars free of internal dissipation. 
Then we will include the pumping and internal dissipation 
only. This will lead to an upper bound on h for gravitational 
radiation reaction free stars. Real stars will be acted upon 
by both gravitational radiation reaction and internal dissi- 
pation. Therefore, at a given frequency, the upper bound on 
h for real stars will be less than the minimum of the two 
separate bounds. In fact, we have already seen that in prac- 
tice internal dissipation is more effective than gravitational 
radiation reaction in damping the wobble (section 2.3). We 
will nevertheless include gravitational radiation reaction in 
this section, partly to quantify just how much weaker it is 
than internal damping, and partly to be able to give an up- 
per bound on the gravitational wave amplitude that would 
apply even if the estimates of internal damping strength of 
section 2.3 are too large by many orders of magnitude. 

We begin by neglecting the internal dissipation. From 
equation (klf) we see that evolves exponentially, increasing 



when 



and decreasing when the inequality is 



reversed. The wobble angle will then either increase to its 
maximum value or decrease to zero: 



= 



for 
for 



Substitution of ( [i^ ) into 
the wave amplitude: 



5G* /crust 

2? r,.e 



1/2 



ump 

(42) 

To.e < T-p ump • (43) 
then gives an estimate of 

(44) 



This is valid for Tpump < Tg^g and so we obtain an upper 
bound on h if we set these two timescales equal: 



/.(Tpump — T(^,f^) — 



5G -^crust 



2c3 



1/2 



(45) 



Now consider the case where the gravitational radiation 
is neglected. Then equations of the form of ( ^ ) and ( |43[ ) 
apply again, but with Tg^g replaced by Td. Substitution into 
(p5|) then gives the estimate 



2nG Ic rust ; ^max 

n = — T (fn • 

c* r Td 



(46) 



This is valid when Tpump < Td and so an upper bound is 
obtained when the timescales are set equal: 



, / N 2nG /crust ; 

/.(Tpump = Td) — T (pn- 



(47) 



4.3 Linear pumping 

Now consider the case of a torque that is fixed with respect 
to the reference plane. Such a torque would lead to a linear 
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variation in the wobble angle, i.e. linear pumping. In this 
case 



1 



1 



(48) 



We will break this up into two separate problems as above. 

When the internal dissipation is neglected we see that 
there is an equilibrium precession angle at which the effects 
of the two torques balance: 



"^pump 

Substitution into equation 
timate for the wave field: 



h : 



5G 

n 'I ^ crust 

2c* 



^1/2 



(49) 

then gives the following es- 



(50) 



Strictly this is an upper bound, as when Tpump is decreased 
below Tgfi the precession angle as given in ( jiol ) becomes 
unphysically large. More precisely, the above equation rep- 
resents an estimate of h for systems with Tg^e/rpump < ^max 
but an upper bound for systems where the inequality is re- 
versed. 

The analogous calculation for radiation reaction free 
systems gives 



leading to the upper bound 

-^crust ' 

Tpump = Tdj = 7 q 



(51) 



(52) 



Clearly the exponential and linear pumping mechanisms 
are very different. The exponential mechanism leads to pre- 
cession only when the pumping torques overcome the radia- 
tion reaction and internal dissipation effects, i.e. either gives 
wobble angles limited by Smax or gives no precession at all. 
On the other hand the linear pumping mechanism will al- 
ways give a finite precession angle, as indicated in equations 
( p9| ) and (|5^). Also, exponential pumping mechanisms can- 
not lead to precession in systems which are originally not 
precessing, while linear pumping mechanisms can. However, 
the upper bounds on h obtained via the two types of pump- 
ing are similar. Comparing (^ ) with (|5^) and (^) with ( |52| ) 
we see that the upper bounds differ only by a factor ^max 
which is of order to unity for low frequencies. 

4.4 Impulsive pumping 

We will describe any process that acts to increase 9 on 
timescales much less than a precession period as an impul- 
sive pumping mechanism. For an isolated star the angular 
momentum during this interval will be constant (apart from 
any negligible amount carried away by a radiation field, un- 
less the star ejects part of its mass!). It therefore follows 
that the kinetic energy of the star must change. On the 
other hand, if the impulse is brought about via an inter- 
action with another body it is possible for only the kinetic 
energy to change, or only the angular momentum, or both. 
Glitches in young pulsars are of interest as candidates for 
isolated impulsive pumping (section]^, while near-body en- 
counters in dense environments are of interest as candidates 
for non-isolated impulsive pumping (section ^). 



5 ACCRETION TORQUES 

Accretion torques are an obvious place to start when looking 
for mechanisms to pump precession. Not only are they capa- 
ble of exciting wobble, but they can also maintain the spin 
frequency of the system, leading to the possibility of long- 
lived constant wave amplitude sources. Indeed, accretion has 
already been investigated as a means of poweri ng gravita- 
tional wave emission via CFS-type instabilities (Papaloizou 



& Pringle 1975 ; Wagoner 1984) and via quadrupole moment 



asymmetries connected with crustal composition variations 
(Bildsten 1998). 

The torque on the central star is the sum of two parts. 
The first is simply the material torque, i.e. due to the accre- 
tion of angular momentum from matter detaching from the 
disk and falling onto the star. The second torque is due to 
the coupling of the star's magnetic field with the disk. The 
net effect of these torques is to spin-up slowly rotating stars 
but impose a maximum spin frequency for fast stars where 
the accretion torque vanishes. It is useful to construct an 
order-of-magnitude accretion spin-up timescale by consider- 
ing the material torque at the magnetosphere radius: 



1.34 X 10* yr 



/ 
IHz 



B 



109 G 



/ M 



-2/7 



(53) 



where M and Me denote the actu al and Eddington accretion 
rates (Shapiro & Teukolsky 1983). Note, however, that this 
will only be equal to the spin-up timescale for slowly rotating 
stars. The accretion torque will vanish when the corotation 
radius lies just outside the inner disk edge. This limits the 
spin frequency of fast stars: 



526 H. i^V'i^] 



3/7 



(54) 



This can be combined with equation (|5^) to give a lower 
bound on ; for a given frequency and accretion rate: 



r„.^> 1.7x lO^yr 



VlHzj \Me)' 



(55) 



If we assume that this torque has a significant compo- 
nent perpendicular to the star's angular momentum vector, 
then the corresponding wobble pumping timescale must be 
/crust// times this (see equation p^ : 

4/3 



> 1.7 X 10^ yr 



/ 
IHz 



M 
Me 



(56) 



We will make use of this when estimating maximum wave 
amplitudes below. 

Almost all analyses of accretion torques have assumed 
axisymmetry, giving rise to either purely spin-up or spin- 
down torques. However, precession pumping requires a 
torque with a component orthogonal to the spin axis. It 
follows that if we wish to find accretion torques capable of 
pumping precession we need to identify situations in which 
the torque itself would be non-constant. Fortunately, we 
would expect this to be the case when magnetospheric ef- 
fects are taken into account. As pointed out by Lamb et 
al. (1975), the accretion rate and therefore also the torque 
depend upon the balance of gravitational, centrifugal and 
magnetic forces. The torque is therefore a function of the 
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plasma angular velocity, stellar angular velocity and stel- 
lar magnetic moment vectors. However, if the star is itself 
precessing the relative orientation of these vectors will be 
modulated. This in turn must lead to a modulation of the 
accretion torque locked in phase with the precession. It is 
precisely this sort of modulation that we would expect to 
lead to a secular evolution in the wobble angle. 

Lamb et al. (1975) demonstrated that such a modula- 
tion can be effective in exciting large amplitude precession 
of a rigid body, with an excitation timescale of the order of 
the spin-up timescale, i.e. of the order of ^ in equation 
(p3|). For our more realistic stellar model this would corre- 
spond to an excitation on a timescale /crust / 1 times shorter. 
A full description of the conditions under which such pump- 
ing can occur would require a detailed understanding of the 
accretion process. Such a description is still not available. 

To re-enforce the potential complexity of this problem, 
the time variation of each vector involved in this problem is 
summarised in table 0. Note that the vector ridisk describ- 
ing the plane of the disk is not well-defined. Far from the 
central star the disk is likely to be planar. However, Vietri 
& Stella (1998) have argued that a non-aligned dipole will 
tend to lift plasma out of the disk plane. In the inner disk, 
viscosity will not be effective in preventing this out-of-plane 
motion. The inner disk will then undergo a forced precession 
with its normal moving in a cone about the star's angular 
momentum. 

A complete description of how this modulates the 
torque on the star has not been attempted previously, and 
lies beyond the scope of this paper. Nevertheless, it seems 
likely that at least part of the torque will be locked in phase 
with the disk precession. This torque would be most effec- 
tive in exciting stellar free precession if the inner disk ro- 
tated at a rate close to the stellar spin frequency, or — more 
accurately — at the frequency (p. However, for all sensible pa- 
rameters the Lense-Thirring precession frequency of the in- 
ner disk is much less than the star's spin frequency. This 
will lead to a rapid averaging of the effects of the torque on 
the wobble angle. It therefore seems unlikely that the forced 
precession of the inner disk plays a crucial role in wobble 
pumping. We will simply note that the notion of a single 
fixed disk orientation is flawed. The vector ndisk in table ^ 
should be regarded as an average orientation, which will de- 
scribe well the outer part of the disk, but not necessarily the 
inner part. Given that the accretion torque may be sensitive 
to the relative orientation of each possible pair of vectors in 
table 1^ the full behaviour is complex, with the torque being 
modulated on a number of timescales. 

To gain a little more insight consider the toy model 
where the torque on the star always points along ndisk, but 
has a magnitude that monotonically increases with the an- 
gle between ridisk and m. Then by considering the relative 
orientation of the vectors of table (|^) it is straightforward 
to find a precessional phase at which the torque, when aver- 
aged over several spin periods, pumps the precession. This 
occurs even when the wobble angle is initially zero. How- 
ever, half a precession period (i.e. n/^) later the relative 
positions of m,nd,J have changed. Then, again when av- 
eraged over several spin periods, it is found that the wobble 
is either damped (if 9 < x) or pumped (when 6 > x)y but 
this pumping is not as strong as before. This is an example 
of an oscillatory torque. However, the monotonicity guaran- 



tees that the pumping over one half of a precessional phase 
exceeds the damping over the next, giving rise to a secular 
increase in 6 which depends on the values of 9 and X- In 
the limit of small wobble angles this increase is proportional 
to 9, giving exponential pumping. Thus, even this simple 
model reproduces a number of the features identified in the 
general arguments above. 

A more sophisticated model is provided by Wang & 
Robnik (1982). Here, the torque on the star is modelled using 
a magnetospheric interaction. The field outside of the star 
is the sum of the usual dipole field plus an additional field 
due to currents induced in the disk. This additional field is 
supposed to be toroidal with respect to the disk's symmetry 
axis and to be significant only near the magnetic poles of 
the star. Wang & Robnik then calculated the torque due to 
the interaction between this toroidal field and the currents 
internal to the star. They found that the torque on the star 
was qualitatively of the same form as the torque that will 
be considered in section ^ when we look at electromagnetic 
torques on isolated pulsars. As will be shown in detail, such 
a torque leads to oscillatory pumping, and also to a secular 
evolution in 9, which for small wobble angles is exponential 
in nature. Again, a simple model has reproduced a number 
of the features identified in the general arguments. 

Of course, the rather general frequency counting argu- 
ments above tell us little about the magnitude of the modu- 
lation. Neither do they say whether the precession is pumped 
or damped. Resolution of such issues requires a detailed 
model of how the torque depends on the angles involved. 
Nevertheless, they show that pumping mechanism could ap- 
ply in accreting systems. We will now consider the particular 
cases of exponential, linear and impulsive pumping and in- 
vestigate the gravitational wave fields they would produce. 



5.1 Exponential accretion torques 

In this case equation ( pl| ) applies, with Tpump = Ta,e. As 
described in section 4.2, we will divide our analysis into two 
parts — one where internal dissipation is ignored, and the 



other where gravitational radiation reaction is ignored. 

Begin by neglecting internal dissipation. Then argu- 
ments identical to those of section (4.2), save for the re- 



placement of Tpump with Ta,e, lead to the upper bounds 



h{T^ = Tg^e) = 1.6 X 10" 



for f < fe and 



/i(-ra 



Tg,0) 



7.4 X 10" 



f 100 Hz 

/ Mbrcak \ / 

100 Hz 



2/3 



1/2 



(57) 



M 
Me 



Ikpc^ 



f 



r 

8/3 



M 



1/2 



(58) 



for / > fe. 

Now neglect gravitational radiation reaction. This leads 
to the upper bounds 



h{T^ = Td) = 4.3 X 10" 



(^) 
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Table 1. This table summarises the temporal behaviour of the vectors of importance in the accretion problem. 
The final column gives their rotation rate as viewed from the inertial frame. As discussed in the text, the disk, 
particularly in its inner regions, may not be planar and stationary, so that n^isk below should be regarded as an 
average disk orientation. Also note that the motion of m is rather complicated — it moves on a cone of half-angle x 
about Tij, while itself moves on a cone of half-angle 8 about J. These two rotations combine to give the motion 
of m indicated. In this case the angular velocities are average values. 



Vector Description Motion Frequency 

'T'disk Normal to disk None 

J Angular momentum of star None 

n Angular velocity of star Cone of half-angle 3edS/2 about J fp 

Tid Axis of star's deformation Cone of half-angle 9 about J </> 

m Dipole axis Cone, |x ~ ^1 < half-angle < |x + ^li </> + '</' for < Xi 

about J (j) for S > X 




log(t) 

Figure 4. The maximum wave amplitude for an accreting star 
when secular accretion torques act. We have assumed «broak = 
10~^, r = 1 kpc, and accretion at the Eddington rate. The curve 
?i(Ta = Tg g) is the bound from balancing gravitational and accre- 
tion torques. The curve /i(Ta = rj) is the bound from balancing 
internal dissipation and accretion torques, with n = 10^. 




for / > fe. 

These curves are plotted in figure m with the values of 
Wbreak, r, M and n indicated by the above parameterisations. 
As can be seen, at all frequencies the signal is limited by the 
internal dissipation, despite the high value of n assumed 
(three orders of magnitude larger than the largest value es- 
timated by Alpar & Sauls 1988). We therefore see that the 
wave amplitude due to exponential pumping by accretion 
torques is almost certainly limited to undetectable values 
by internal dissipation. We will not consider this mechanism 
any further. 



5.2 Linear accretion torques 

In this case equation ( ji^ ) applies, with Tpump = Ta,e- This 
will lead to upper bounds on h that coincide with the low 
frequency portions of the curves in figure ^. Therefore in- 
ternal dissipation still limits the amplitude to undetectable 
values. 



5.3 Impulsive accretion torques 

We will now discuss a pumping mechanism that will ap- 
ply even in the absence of the gating described previously. 
Suppose the accreting plasma is clumped, and is not en- 
tirely confined to the plane perpendicular to the star's ro- 
tation axis. Various instabilities are likely to lead to such 
a situation, as discussed in Lamb et al. (1985). Indeed, the 
clumping of accreting plasma is a key ingredient of the 'beat- 
frequency' models. These ascribe the quasi-periodic oscilla- 
tions observed in many X-ray systems to a variability in 
the accretion rate. This variability is due to a beating be- 
tween the frequency of the innermost stable orbit and the 
star's spin frequency. It therefore seems likely that the phe- 
nomenon of clumping in accreting systems is generic. 

Denote by Ati the (average) timescale in which a clump 
transfers its angular momentum to the star. The wobble 
angle will change impulsively due to the accretion of such 
clumps providing the transfer of angular momentum takes 
place over a narrow range of rotational phase, i.e. provid- 
ing Ati ^ P, where P denotes the star's spin period. If 
in addition the clumps arrive at random time intervals of 
average At2 then, providing Pfp <^ At2, a random walk in 
6 will occur, even though no gating occurs. In the absence 
of gravitational radiation reaction and internal damping the 
wobble angle will then grow as t^^^. This is clearly a less ef- 
fective pumping mechanism than the linear one considered 
previously, and so will not lead to interesting gravitational 
wave amplitude. 



6 ELECTROMAGNETIC TORQUES 

Models of the electromagnetic torque on a spinning neu- 
tron star have split into two classes — the Goldreich & Julian 
(1969) model where a dipolar magnetic field embedded in a 
the star generates a strong electric field at the stellar surface 
which rips out charged particles. These then propagate along 
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the magnetic field lines, forming a magnetosphere. This ap- 
plies even when the dipole and rotation axes coincide. The 
radiation the charged particles emit then carries energy away 
from the star, implying the existence of a braking torque. 

The second class does not have a magnetosphere, and 
instead models the star as a perfect conductor with a mag- 
netic dipole embedded, surrounded by vacuum (Ostriker & 
Gunn ] |969[ ). If the dipole is inclined to the rotation axis, 
electromagnetic radiation at the spin frequency is emitted. 
This provides the braking torque. Such models are clearly 
oversimplified. However, the magnetosphere models are so 
complicated that very littl e progress ha s been made in un- 
derstanding their details (Michel 1991); questions such as 



how the electromagnetic torque depends upon the angle 
between the spin axis and dipole axis remain unanswered. 
Therefore, in this section we will look at a magnetosphere- 
free model in which the torque has been calculated as a func- 
tion of the spin and dipole moment vectors. Even though the 
details of such a model may prove to be incorrect, it is likely 
that it will reproduce qualitatively many features of real 
stars, and in particular identify the important timescales 
involved. 

Most discussions of magnetosphere-free torques are 
based on the model of Deutsch (1955), who considered a 
perfectly conducting sharply bounded sphere, rigidly rotat- 
ing in a vacuum. The internal magnetic field was assumed 
to be symmetric about some axis inclined to the axis of ro- 
tation. The external magnetic and electric fields were then 
calculated. These fields were then used by Davis & Gold- 
stein (1970) and Michel & Goldwire (1970) to calculate the 
torque on the star. These authors investigated the effect of 
this torque on spherical stars. They were able to show that 
the torque caused the magnetic dipole to align with the ro- 
tation axis on the electromagnetic spin-down timescale. 

Goldreich (1970) extended the analysis to examine the 
effect of torques on precessing bodies. He showed that when 
the free precession period was less than the electromagnetic 
spin-down timescale the non-sphericity completely altered 
the body's evolution — the magnetic moment no longer aligns 
with the spin axis. Instead the torque serves to either pump 
or damp any free precession the body may be undergoing, 
depending upon whether the angle x between the dipole mo- 
ment and the body's deformation axis is greater or smaller 
than sin~^(2/3)^''^ ~ 55°, respectively. We will refer to this 
process as the Goldreich mechanism. This is very different 
from the gravitational radiation re action case, in which the 
wobble motion is always damped (Cutler & Jones 200C). 

The issue of calculating the electromagnetic torque was 
taken up again by Good & Ng (1985). Much of their analysis 
was concerned with Goldreich- Julian type magnetosphere 
torques. However, in the course of their calculation they 
found an error in the Deutsch fields. When this is taken 
into account the Deutsch vacuum torque T is given by: 



2Q 1 



(n-m)(nxm)(61) 



where fl denotes the angular velocity, R the stellar radius 
and m the dipole moment. (The only difference between 
this and the torque as calculated from the original Deutsch 
paper is a factor of —1/5 in T2). 

The first term, Ti, has components both parallel and 
perpendicular to the spin axis, and scales as Q,'^. It is known 



as a non-anomalous torque. It causes the energy and angular 
momentum of the star to decrease. Its component parallel to 
J is responsible for spin-down, while its component perpen- 
dicular to J will bring about a secular change in the wobble 
angle. 

The second term, T2, is exactly perpendicular to the 
angular velocity and scales as . It is known as an anoma- 
lous torque, and is caused by the near-zone fields. It does 
not lead to a loss in energy or angular momentum. This 
torque does affect the wobble angle and spin rate of a freely 
precessing star, but its effects average to zero over one free 
precession period. 

Following equation (^l|) we will define two timescales to 
characterise the torque. The first corresponds to Ti and is 
simply the spin-down timescale f /f — ZIc? /2E!' (jp' B^: 



To = 3.64 X 10"^ yr 



10^" G 
B 



^ ^lOOHz^ ^ 



f 



sin^ X 



(62) 



The corresponding timescale on which this torque would 
evolve the wobble angle is /crust shorter than this. In the 
above formula, and those that follow, we will place the geo- 
metric factors that follow from the particular form of vector 
and scalar products of equation ( |6l| ) in square brackets. In 
this way it will be clear what conclusions can be drawn if 
only the approximate timescales of the model hold, i.e. if in 
more realistic magnetosphere models the above geometric 
factors were found to be incorrect. Such a situation would 
correspond to replacing the term in square brackets with 
some other geometric factor. 

The second timescale corresponds to T2 and we define 

as 



76.2 yr 



10^^ G 
B 



100 Hz 



sin^ X 



(63) 



It is a factor ^IR/c shorter than Tc. Again, the correspond- 
ing timescale on which this torque would evolve the wobble 
angle is /crust// shorter than this. 

Thus the effects of the two torque terms are qualita- 
tively different. The term Ti results in a secular variation 
in the free precession angle and spin rate on the electro- 
magnetic spin-down timescale. The term T2 produces no 
such secular variation. Instead it causes oscillations in the 
spin frequency and wobble angle on the (much shorter) free 
precession timescale. We will refer to this as a non-secular 
mechanism. We will examine both torques in terms of their 
effect on gravitational wave generation below. 



6.1 Non-secular electromagnetic torques 

The torque T2 has received very little attention previously, 
undoubtedly because it does not cause any secular varia- 
tion in the star's motion. However, as can be deduced from 
equation ( |6l| ) its magnitude and orientation with respect to 
the reference plane varies over one precession period, as m 
rotates with respect to this plane. Recently Melatos (1999) 
made use of this torque to model spin-down irregularities in 
magnetars, modelled as rigid bodies. He considered the case 
where the timescale fc was of similar duration to the free pre- 
cession period: fo ~ Pfp. This similarity of timescales gives 
a rather irregular or 'bumpy' spin-down rate, which Melatos 
calculated numerically. We will take a somewhat simpler 



© 0000 RAS, MNRAS 000, 000-000 



12 D. I. Jones and N. Andersson 



view — when fe is longer than Pfp the effect of the torque 
will be to simply cause a smooth variation of free preces- 
sion parameters (e.g. 9 and (p) calculable using perturbation 
theory. Even when the two timescales become comparable 
we would expect our results to apply to order-of-magnitude 
accuracy. 

Substitution of T2 into (|l^) shows that for small wobble 
angles, 9 varies sinusoidally, giving oscillatory pumping of 
the form discussed previously. An estimate of the magnitude 
of this angle follows from the timescale fe described above: 

9ns ~ ^7^^. (64) 

'^c J crust / J 

From equation (^o|) this leads to a wave field 



h 



Gulp 
2c* r ' 
This gives 



/i = 1.8 X 10" 



{"¥) m 



Stt tan X 



(65) 



(66) 



where we have parameterised in terms of Tc rather than fe. 
As above the geometric factor deduced from the full pertur- 
bative calculation has been separated from the rest of the 
formula. 

Clearly this is extremely weak. The amplitude increases 
as the electromagnetic spin-down timescale decreases, so we 
should focus attention on fast spinning strongly magnetised 
stars, i.e stars very soon after birth in supernovae. Such stars 
will be hot, with temperatures of order lO'^^ K immediately 
after birth. However, this mechanism can become active only 
when the outer phase has solidified to form a crust, and so 
we can apply equation (|6^) only to stars which cool to the 
melting temperature in a time less than Tc. 

As described in Haensel (1997), the crust will not have 
a single well defined melting temperature, as the deeper 
parts (with density ~ 10^'^gcm^) will melt at tempera- 
tures just under 10^" K, while the outer crust (with density 
~ 10^^ gcm^) will melt at temperatures just above 10^ K. 
There is therefore some ambiguity in what to take as an 
average melting temperature, with a corresponding ambi- 
guity in the timescale for crust formation. For the range of 
temperature identified, an upper bound woul d be one year, 
with a lower bound much shorter than this (Haensel 1997 
Shapiro & Teukolsky 19831). 



If the solidification timescale is indeed of order one year, 
equation (^) then leads to the tiny signal of order 10~^^ 
for a star born at the Galactic centre. To be detectable over 
a one-year integration the geometric factor would have to 
exceed lO'^, corresponding to very small dipole inclination 
angles, x < 2°. 

If the star's outer layers were to cool more rapidly than 
this we can consider the effective amplitude found by mul- 
tiplying the wave amplitude by the square root of the num- 
ber of revolutions performed in one electromagnetic braking 
timescale: 

■lyrNi/2 



/..ff = 2.5xl0-«(i^) (if) 

(1 



/ 





'3 1 " 




5n tan^ x 



.(67) 



,100 Hz. 

Even with fast (less than one-year) cooling, a star born at 



the Galactic centre remains undetectable unless the geomet- 
ric factor amplifies the signal. A dipole inclination of 5° gives 
an amplification of 25. At 100 Hz the signal is detectable for 
Tc — 1 month, and a polar magnetic field of 2 x 10^^ G. 
Such an event could correspond to the birth of a magne- 
tar. However, the event rate for such an occurrence is very 
low — surely less than one a century, so that the probability 
of such a source being born during the operational lifetime 
of Advanced LIGO is small. 

The torque T2 is due to the near-zone fields and does 
not lead to energy being radiated to infinity. Therefore it 
is conceivable that it may continue to act in an accreting 
system. In this case the spin frequency would be maintained 
by the accretion torque. To obtain an upper bound on the 
wave field at a given frequencyset To equal to its minimum 
value as given by combining (M) and (p2): 



: 3.4 X 10^ 



VlOO Hz/ 



M 



(68) 



Substituting this into equation (^6|) gives a maximum wave 
amplitude of 



Vax(/) = 2.4 X W-'' (i^) 



V / 



lkpc\ / 100 Hz V''^ / M 



Me 



(69) 



This is tiny even in comparison to the wave amplitudes from 
isolated stars. This is because in accreting systems the two 
factors which give large wave amplitudes — high spin fre- 
quency and high polar magnetic field strength — compete. 
Highly magnetised stars have large Alfven radii and there- 
fore low equilibrium spin rates. 

To sum up, the non-secular electromagnetic torque does 
not seem to be a good pumping mechanism for gravitational 
wave generation. Isolated stars at the Galactic centre with 
high spin frequencies and large magnetic field strengths can 
produce detectable signals if the geometric factors which en- 
ter the calculation are favourable. However, if born in super- 
novae, their rapid spin-down rates would require extremely 
fast (less than one-year) cooling of the outer phases of the 
star. Accreting stars where this near-zone torque continues 
to act produce even smaller gravitational wave signals, as 
the requirements of high spin frequency and high magnetic 
field strength oppose. 

6.2 Secular electromagnetic torques 

Having investigated the non-secular effect of electromagnetic 
torques on gravitational wave generation we will now con- 
sider the possibility that such torques may pump precession 
in a secular way. Goldreich (1970) demonstrated that the 
torque Ti of equation (^l|) does indeed lead to such evolu- 
tion for a rigid body. This generalises at once to our elastic 
shell/fluid core model. Substitution of Ti into equations ( |lll ) 
and leads to a set of two coupled differential equations 
connecting and 9. Providing the free precession period is 
less than the spin-down timescale we can average over a free 
precession period to give: 



-Q(/)^ sin^ X 



-^crust 

where 



3 . 2 



(70) 
(71) 
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2w/ 
3/c3" 



(72) 



(Small terms proportional to A/d//crust and 6 have been 
neglected). Then the spin-down follows the usual power law: 

-1/2 

(73) 



1 



2 sin^ xi 



where Tc is defined in equation (^2|). The wobble angle 
evolves according to 

A 



2 sin^ xt 



where 



A = 



I sin^ X 



2 sin'' X 



(74) 



(75) 



As can be seen, 6 increases or decreases depending upon the 
sign of I sin'^ X ~ li s^nd increases most rapidly for x = 90°, 
in which case 



0{X = 90° 



(76) 



In reality gravitational radiation reaction and internal 
dissipation will act on the spinning down star also, and so we 
should include their effects. Then equations ( [to] ) and (^]) ac- 
quire extra terms. The gravitational spin-down of equation 
is added to ([70[), while the gravitational alignment of equa- 
tion (|l^ is added to ([ti]). The internal dissipation alignment 
rate corresponding to equation (|l^ is added to ([ti)). 

We will present the x = 90° results where the electro- 
magnetic pumping is greatest: 



-ad) 



2 Icvus 

where 



/crust 27rn 



(77) 
(78) 



(79) 



The terms proportional to a are due to electromagnetic 
torques, those proportional to /3 to gravitational torques, 
and the term proportional to n^^ describe internal damp- 
ing. In terms of timescales equation (Fsl) would be written 



1 



1 



1 



(80) 



where for convenience we have introduced the electromag- 



netic whhble pumping timescale 



re,e 



= 2re 



(81) 



Note that the timescales in this problem are functions of 
the decreasing frequency, so this equation will generally not 
have simple exponential solutions. 

We will apply this formalism to a young rapidly spin- 
ning down neutron star. Its initial wobble may be due to 
glitches arising from the rapid spin-down (section |^, or due 
to some other process. We will not concern ourselves with 
this issue here. We will simply examine the evolution of the 
precession angle assuming that the Goldreich mechanism 
acts. 



Fortunately it is possible to simplify the above equa- 
tions somewhat. Firstly the gravitational radiation reaction 
timescale for spin-down is much longer than for gravitational 
alignment, so in any situation where the timescales in ( |80[ ) 
are similar the gravitational spin-down term is negligible, i.e. 
the last term of (|7^) may be neglected. Also, the damping 
due to gravitational radiation reaction will be much weaker 
than the damping due to internal dissipation. Then the sec- 
ond term of (|7^) may be neglected. 

For such a star the evolution in <f> is driven only by 
the electromagnetic torque and equation ( [t^ ) applies. The 
evolution in 6 is more complex. Integration of ([TSf) depends 
on the variation of Aid with the frequency. We will consider 
the case of Coulomb deformations, and assume Aid oc <f>^ 
(i.e. that the crust is relaxed). We then find 



1 + ' 



exp 



6rd 



1 + ^ 



-1/2 



- 1 



(82) 



We can now build up a picture of the evolution of a star 
which is acted upon by a secular electromagnetic torque of 
this form. The scaling of the timescales with frequency are 
crucial: Td oc l/AId4> oc l/4>'^, while Tc oc l/<p. If the star 
is born spinning sufficiently fast the steeper dependence of 
Td on frequency will give Td ^ Ta and the above equation 
reduces to 

t 

Td - 

As the star continues to spin-down the two timescales be- 
come comparable an d ev entually 9 will begin to increase. At 
late times equation (S.2) reduces to 



9 K 9o exp 



(83) 



) exp 



(84) 



This is of the same form as ([76|), save for the extra expo- 
nential factor. This factor can be interpreted as the extent 
to which the wobble was damped early in the star's life. An 
example of this behaviour is shown in figure ^. Time is mea- 
sured in units of Tc and we have chosen Td = O.IStc. An 
example of such a star where Coulomb forces provide the 
distortion could have b = lO"'', /crust// = 0.1, B = lO" 
G and an initial spin frequency of 200 Hz. Then To = 9 
years. The initial exponential decrease in 9, on a timescale 
of approximately one year, is clearly seen. Of course, this ne- 
glects the possibility that a CFS-type instability (such as an 
r-mode) might ra pidly break the star's rotation (Andersson 



Kokkotas 2001) 



In this way electromagnetic radiation reaction could 
provide a way of setting young neutron stars into free pre- 
cession. Given some initial non-zero wobble angle and the 
condition re,9 < Td, the wobble angle increases. Of course, 
this increase will end when 9 reaches the maximum value 
as set by the crustal breaking stress. From that point on 
the wobble angle will remain fixed at its maximum value, 
while the spin frequency continues to decrease according 
to equation ([ts}). The gravitational wave amplitude is pro- 
portional to AId4>^ . Given our assumption A/a oc , this 
leads to a wave amplitude proportional to 0*, which in turn 
evolves as (l-|-2f/ro)~^. In other words, for all but the earli- 
est stages of evolution the gravitational wave amplitude de- 
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Figure 5. The evolution in the wobble angle 9 under the ac- 
tion of combined internal dissipation and electromagnetic torques. 
We have put re = 0.13 X r^ g. At early times the internal 
dissipation dominates and 9 decreases exponentially. At later 
times the electromagnetic torque dominates and 9 increases as 
(l + 2t/rc)^/«"-'. 



creases steadily with time — we have been able to maintain 
a non-zero wobble angle only at the expense of introducing 
a powerful spin-down torque. The gravitational wave ampli- 
tude decreases on a timescale of order Te. In section |^ we 
will look at the wave field due to a population of young un- 
magnetised isolated stars, set into free precession at (or very 
soon after) birth, whose precession angles then decay on the 
internal damping timescale ra- Given that the wave field of 
the magnetised stars considered in this section decays on the 
timescale Ts,e < Td, the magnetised stars are not as easily 
detectable, and so we will not pursue them further. Rather, 
in section ^ we will consider in detail the signal strength 
and population statistics of young unmagnetised stars, and 
find that for realistic values of the damping parameter n, 
detection of such a population is not likely. 

To summarise, this non-secular electromagnetic pump- 
ing can increase the wobble angle of an isolated star. How- 
ever, once the wobble angle has increased to its maximum 
value the associated spin-down torques dominate the evolu- 
tion in the wave signal, with the result that the amplitude 
decreases on the electromagnetic spin-down timescale. 

This motivates us to consider the possibility that secular 
electromagnetic pumping remains active even in an accret- 
ing system: The accretion torque could then maintain the 
star's spin frequency, while the electromagnetic torque main- 
tains the free precession. It is more difficult to justify this 
combining of accretion and electromagnetic torques than in 
the case of the non-secular electromagnetic torque, as a sec- 
ular variation in 6 changes the energy and angular momen- 
tum of the star, and so must correspond to fiuxes of energy 
and angular momentum to infinity. Certainly it is clear that 
the exact form of the torque Ti will no longer apply, as its 
spin-down part at least will be suppressed by the accretion 
environment. We will therefore dispense with the details of 
the above, i.e. the x ^^^^ numerical factors, and will assume 
that the star is in spin equilibrium as described by standard 
accretion theory. Then the star is modelled by 



= 0. 



(85) 



(86) 



The timescales in ( pq ) are constant by virtue of (pq), 
and so the solution will be exponential in form. However, 
we now have to include the details of the accretion spin-up 
torque, the electromagnetic pumping torque, and internal 
dissipation. In fact, it is possible to combine formulae de- 
scribing these three phenomena to obtain an upper bound 
on the gravitational wave signal independent of the mech- 
anism producing the deformation. Firstly, we can combine 
equation ([64 ) for the electromagnetic spin-down time with 
equation (M) which gives the maximum spin frequency of a 
star with a given magnetic fields strength and a given ac- 
cretion rate. This provides the lower bound on Tc given by 
equation (^). However, for the Goldreich mechanism to be 
operative we must have Te,e ~ 2rc /crust// < Td- The above 
equation can then be used to give a lower bound on Td, at 
a given spin frequency, accretion rate and value of /crust//. 
But Td oc /crust /A/d, and so we obtain an upper bound 
on A/d. This immediately become an upper bound on the 
gravitational wave amplitude. Carrying out the arithmetic 
we find: 



h{r^ = Td) = 5.0 X 10" 



( Mf J ( 



/ 



Mb / V 100 Hz 



2/3 



for / < fo and 
h{Ta = Td) = 2.4 X 



,„-31 /Wbrcak\ / _M \ f 100Hz\ 

V 10-3 ) [m^J [ f J 



(87) 



4/3 



( — 

Vio^ 



for f > fe- This is a tiny amplitude. Even for n — 10^ it 
reaches a maximum of only 10"'^" , at a frequency of 69 Hz for 
the values of Ubrcak chosen. Therefore the gravitational wave 
amplitude of accreting stars acted upon by electromagnetic 
torques of this form is limited to uninteresting values by 
internal dissipation. 

To conclude, secular electromagnetic torques are not 
likely to lead to detectable levels of gravitational radiation. 
In the case of isolated stars the torque may be able to in- 
crease an initially small wobble angle to its maximum value, 
but only at the expense of introducing a strong spin-down 
torque. Even if this part of the torque were to remain active 
in accreting stars the wave amplitudes are uninteresting. 
The bound on h obtained by considering the competition 
between the electromagnetic torque and the internal dissi- 
pation timescale leads to a bound on h much lower than 
even the Advanced LIGO sensitivity. We will therefore not 
consider the effects of electromagnetic torques any further. 



7 NATAL PRECESSION AND GLITCHES 

Given their violent birth in supernovae it is tempting to ex- 
amine the possibility that neutron stars are set into preces- 
sion when born. As discussed in section the high tem- 
peratures generated by the implosion will lead to entirely 
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fluid stars. Only when the outer layers have cooled will a 
solid crust form. It is possible that when the crust is in the 
process of solidifying a substantial excitation of a mode of 
oscillation of the fluid exists. This excitation could simply 
be due to the supernova explosion itself, or could be due to a 
CFS-type instability. Such a scenario has already been inves- 
tigated in the context of r-modes, where the crust formation 



-i9r 



was a hindrance to gravitational wave emission (Lindblom 



Owen & : Ushomirsky 200C ; Andersson et al. 200C ) . The exact 



outcome is not clear, particularly as the crust will not so- 
lidify at all points simultaneously, but will form flrst where 
the fluid velocity is smallest. Given the complexity of the 
process, the possibility of the crust/fluid core system being 
created in a precessional state cannot be ruled out. 

Even if the star has settled down into an axisymmetric 
configuration at the time of crust solidification, it is possible 
that it might be set into free precession soon after. As noted 
by Lyne (1996), the phenomenon of glitching — the sudden 
increase in spin frequency of a pulsar — is most common in 
young pulsars. This phenomenon certainly requires a solid 
crust. Therefore, a young neutron star that has cooled suffi- 
ciently will begin to glitch. Also, some theories of glitching 
associate the sudden change in spin frequency with a frac- 
ture in the crustal structure (see Ruderman (1976, 1991a,b) 
and Link, Franco & Epstein (1998) for details). If this change 
results in a shift in the principal axis of the moment of in- 
ertia tensor the star will precess ( Pines fc Shaham 1973a ; 
' Pines fc Shaham 1973b|; [Link, Franco fc Epstein 1998|). In 



this way it is possible that stars may acquire a precessional 
motion very soon after birth. 

We therefore wish to investigate the gravitational wave 
background due to a population of young spinning-down 
neutron stars that were set into precession at, or soon after, 
birth. We will describe this as natal precession. Although we 
have identified fluid modes and glitches as possible ways of 
producing this natal precession, the following analysis would 
apply regardless of the source of the wobble. 



7.1 Gravitational wave amplitudes 

We will not include any pumping mechanisms in our analy- 
sis, so that the wobble angle simply decays under the com- 
bined influence of internal dissipation and gravitational ra- 
diation reaction. Also, we will begin by considering stars for 
which electromagnetic spin-down torques are negligible. We 
would then have 



1 



1 

Td 



(89) 



If we set the initial wobble angle equal to its maximum value 
we have 



9{t) = Smax exp 



(-7) 



(90) 



for t > and is zero for t < 0. 

In order to gain insight into the detectability of the 
gravitational wave field due to such a source we shall assume 
that a matched filter can accumulate signal only for an inter- 
val T or for an interval of one year, whichever is shorter. Fig- 
ures P and ^ show the effective amplitude for a star deformed 
by crustal Coulomb forces, with a strain angle "Ubroak ~ 10""^, 
at a distance 1 kpc. Each figure plots h for three different 
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Figure 6. The effective amplitude from an isolated star which ini- 
tially precesses at its maximum wobble angle with damping due to 
gravitational radiation reaction only. We have put Wbrcak — 
and r = 1 kpc, and have assumed that the matched filter accu- 
mulates signal for an interval t defined in equation (^) or for an 
interval of one year, whichever is shorter. The deformation is due 
to Coulomb forces, with rigidity parameters b = 10""^, 10"* and 
10~^ as indicated. The knee that appears on the b = 10~^ curve 
is due to the gravitational radiation reaction timescale falling to 
below one year at high rotation rate, limiting the signal accumu- 
lated by a matched filter. 
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Figure 7. The effective amplitude from an isolated star with the 
same parameters and assumptions as in figure (^), except now 
n = 10^. This corresponds to the upper limit on n as estimated 
by Alpar & Sauls (1988). 



values of the rigidity parameter: b = 10~"',10~^ and 10~^. 
We have assumed a reference oblateness equal to that of 
a rotating fluid, i.e. eo = en, so that Aid = 3/star6en/2. 
In figure |6| on ly the gravitational radiation reaction term of 
equation (p9[) is included, i.e. this is the n 00 limit. Figure 
has n = 10'* , which is the upper bound on n based on the 
arguments of Alpar & Sauls (1988). This level of damping 
reduces the signal significantly as compared to figure ^ for 
frequencies in excess of 100 Hz. For example, a b — 1Q~^ 
star is now barely detectably by a first generation interfer- 
ometer, while the more plausible b = 10~^ star is now barely 
detectable by an Advanced LIGO. 
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7.2 Statistical arguments 

Such plots as these are useful as they show how the strength 
of the internal damping affects the strength of the gravita- 
tional wave signal from an isolated star. However, in order 
to decide whether such sources are of interest it is necessary 
to consider the related issues of their event rate and distance 
from the Earth. To do this we need to consider the statistics 
of neutron star births. 

The analogous problem of the gravitational signals due 
to a population of young non-precessing triaxial neutron 
stars spinning down due to gravitational wave emission has 
been considered by Blandford as reported in Thorne (1987). 
Blandford made use of the following argument: Consider a 
source born in a supernova explosion that is spinning down 
on a timescale r. Then there will be r/AfsN such sources 
with age r or younger in the Galaxy, where AtsN is the inter- 
val between Galactic supernovae. Blandford then modelled 
the Galaxy as a flat disk of radius R, which allowed him to 
estimate the distance to the nearest such source. He then 
combined this result with gravitational spin-down and wave 
amplitude formulae from the quadrupole formalism to show 
that the gravitational signal arriving at Earth from the near- 
est such source is proportional to \fKta^ and is independent 
of the source's frequency or triaxial ellipticity. 

We wish to make the analogous argument for a pop- 
ulation of young isolated precessing neutron stars. We will 
follow Blandford and consider stars where electromagnetic 
torques are not important. We have already shown in sec- 
tion 3.2 that electromagnetic torques can, if the Goldreich 



h = 2Ax 10" 



pumping mechanism is active, cause young neutron stars to 
be set into free precession, but only at the cost of intro- 
ducing a powerful spin-down torque, and so a population of 
unmagnetised stars is probably more easily detectable than 
a population of magnetised ones. 

The relevant timescale t on which the gravitational 
wave signal decreases is no longer the gravitational spin- 
down timescale, but is instead the free precession alignment 
timescale, dominated in almost all physically plausible sce- 
narios by internal dissipation. This means that Blandford's 
triaxial result where /i is a function of AfsN only no longer 
holds. 

We will model the Galaxy as having a radius R and 
thickness D, giving a volume of order R?D. This contains 
r/AtsN stars of age r or younger. The average separation 
of this population of young stars is then {AtsT<sR^ D /tY^^ . 
Explicitly 



Ar = 1.4 kpc 



f Ats 



30 yr 



1/3 



10' yr 



1/3 



(91) 



where we have put R — 10 kpc and D — 1 kpc. Of course, 
a more accurate model would take into account the rate 
of star formation as a function of Galactic position, with 
different rates applying in the central bulge and spiral arms, 
for instance. Nevertheless, equation (^) represents a useful 
first approximation. 

Equations (^^ and (|3^) give the wave amplitude of a 
precessing star in terms of r,f and t^. We can then set r 
equal to Ar as given in equation (^l|). The wave amplitude 
thus obtained will, subject to statistical variation, be the 
field at Earth due to the closest source of age ra or less. In 
full: 



10-* V 100 Hz/ 
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for f < fe and 



h=l.lx 10" 



/ ^tbrcak \ 



^ 100 Hz 



/ 



/30yr\ 
VAtsN/ 



1/3 



10' yr 



2/3 



(93) 



for f > fe- We are considering stars with r in excess of 30 
years, so these wave amplitudes can be compared against the 
noisecurves for 1-year matched filter integrations, e.g. the 
noisecurves of figure (^. As can be seen, signals oi h = 10~^^ 
lie above the Advanced LIGO noisecurve for a wide range of 
frequencies (25 400 Hz) which may well include the initial 
spin frequencies of neutron stars. However, we are consider- 
ing a 'blind search' i.e. an all-sky search without any prior 
knowledge of the source's position and little idea of its spin 
frequency. In order to minimise numerically-generated false 
alarms we will therefore take our detection criterion to be 
h > 10"^''. Also we will set Wbrcak ~ 10~'. For definiteness 
consider the case where all stars are born with a spin fre- 
quency / = 100 Hz. If we consider the optimistic case where 
internal dissipation is neglected (n — > oo) then the above 
equations show that the nearest such star is detectable for 
r < 30 years. This is on the edge of the applicability of our 
simple statistical model — there would only be one such star 
in the Galaxy! The probability of detecting such an object 
during the lifetime of an Advanced LIGO detector would be 
borderline. At follows at once that damping parameters of 
less than n — 10^ would be unlikely to lead to a detectable 
population of isolated stars. 

If, when the advanced detectors go on line, advances in 
computer power, detector noisecurves or search algorithm 
design permit a more sensitive search down to ^ = 10"^^, 
then we obtain Td = lO' years for detection. There would 
be a population of around 30 stars of this age or less in 
the Galaxy. Equation ( pl| ) shows that the closest would be 
about 1 kpc from the Earth. Figure ^ shows that in the case 
of Coulomb deformations such a star would have to have 
a rigidity parameter of 10~*, an order of magnitude larger 
than the canonical value. 

We are now in a position to summarise our results. We 
have considered the wave field due to a population of isolated 
stars set into precession with 6 = ^max soon after birth. The 
signal following matched filtering is shown in figures |^ and 
1^ for the case of Coulomb deformations with Ubroak ~ 10~' 
and r = 1 kpc. They suggest that even for damping as strong 
as n = 10*, a star with rigidity parameter b — 10~* is de- 
tectable out to about 1 kpc by an Advanced LIGO. How- 
ever, when the statistics of neutron star births are included 
it is found that with this level of internal damping it is un- 
likely that, at a given time, there exists such a precessing 
star in the Galaxy. However, the situation is very different if 
the internal damping is much weaker. In the limit where it 
negligible, with a detection threshold of 10~^^ it is border- 
line as to whether or not there will exist a detectable star 
in the Galaxy. If this threshold can be decreased to 10"'^'', 
then there could exist a detectable star as close at 1 kpc 
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to Earth. If deformed by Coulomb forces, such a star would 
require a rigidity parameter of 10"* at least. 



8 COLLISIONS IN DENSE ENVIRONMENTS 

The possibility of detecting gravitational waves from neu- 
tron stars set into free precession following collisions with 
other stars has recently been discussed in the literature, de 
Araujo et al. (1994) first pointed out that the high stellar 
densities of globular clusters could lead to a high rate of 
collision between neutron stars and other stars. Given that 
globular clusters contain an excess of millisecond pulsars, 
the latter authors argued that if such collisions were effec- 
tive in exciting precession, interesting levels of gravitational 
wave emission would occur. 

This idea was taken up by Velloso et al. (1997) who 
estimated possible gravitational wave amplitudes at Earth 
based on the millisecond pulsar data for the globular clus- 
ters. However, the formulae employed were taken from de 
Araujo et al. (1994) who effectively assumed that the stars 
were rigid, with an oblateness en. This gave h cc (jP'enO 
and Tg.e oc I/eq. At a given wobble angle and spin fre- 
quency this led to a wave amplitude too large by a factor 
of en/ed ~ 1/6 and an alignment timescale too fast by the 
factor (en/ed)^ ~ l/b^ . 

In this section we will reconsider the issue of collision- 
induced precession. In the terminology of section^ this is an 
example of an impulsive pumping mechanism. This problem 
is very similar to that considered in the last section, natal 
precession, as we are again considering the wave field due 
to stars set into precession at positions and times which 
can only be described in a statistical way, and then spin- 
down and align. Indeed, the problem again breaks down 
neatly into two parts. The first concerns identifying par- 
ticular mechanisms, i.e. types of collision, that lead to an 
interesting level of free precession. The second concerns find- 
ing an event rate for such an occurrence, so that a statisti- 
cal statement can be made regarding the likely detectability. 
Note that as the nearest cluster is of order 1 kpc from Earth, 
equations ( |36[ ) and ( |37|) show that a safe condition for de- 
tectability is Tg^e ^ 10^ years, with the very optimistic as- 
sumptions n = 10^ and a detection threshold ol h — 10"'^'^. 
Additionally, we see from figure ^ and equation (^4|) that 
for a star with h = 10~^ and Ubrcak ~ 10"'' spinning at sev- 
eral 100 Hz, a wobble angle of order 10""' is required for a 
signal-to-noise of 10. We therefore need to identify types of 
collision involving at least one recycled neutron star which 
produce precession angles of order 10"'' and have an event 
rate of order 1/1000 yr~^. These arguments are somewhat 
simpler than those used in the last section, but will suffice 
for this section, as we will find very small event rates for 
collision, so that it will be clear that we can rule out colli- 
sions as a mechanism for generating detectable gravitational 
waves. 



8.1 Collisions of neutron stars with non-compact 
stars 

We begin by considering the collision of a neutron star with a 
non-compact star. Such collisions have been modelled exten- 
sively, e.g. Davis, Benz and Hills (1992), who model neutron 



star-main sequence star collisions, and also neutron star-red 
giant collisions. In the main sequence star case they find 
that a system consisting of a neutron star surrounded by a 
thick accretion disk is formed when the separation at peri- 
astron is ^ 1.75 times the main sequence star's radius. In 
the red giant case they find that a common-envelope system 
is formed when the separation at periastron is ^ 2.5 times 
the red giant radius. For periastron separations significantly 
greater than these values they find that the perturbation of 
the non-compact star is minimal. 

Despite the violent effect such near-body encounters 
have on the non-compact star it is difficult to see how the 
neutron star would be set into fre e pr ecession by such a col- 
lision. As will be shown in section 8.2, the gravitational tidal 
torque on the neutron star due to the non-compact star is 
negligible. This leaves only the material torque on the neu- 
tron star, which will be determined by accretion fiow onto 
its surface. This will be described by the standard theory, 
regardless of the unusual source of the accreting material. 
The accretion rate will be limited to the Eddington value in 
the usual manner, so the torque will not be impulsive. 

It follows that although collisions between neutron stars 
and non-compact stars are important in terms of the popu- 
lation evolution of globular clusters, they are not of use as a 
mechanism for free precession gravitational wave production 
as it is impossible to identify a way in which the collision 
would set the star into free precession. 



8.2 Neutron star-neutron star encounters 

We will now consider encounters between two neutron stars. 
Clearly, if a direct collision were to occur, free precession 
would be the last gravitational wave mechanism that we 
would wish to consider. We will therefore model a near col- 
lision, where both gravitational and magnetic effects will 
come into play, but there is no direct mechanical contact 
between the stars. We will begin by considering the grav- 
itational interaction in a simple Newtonian way. Suppose 
one star has a spin-angular momentum J and centrifugal 
deformation A7n. Then the other star will exert a torque 
on this bulge, causing a forced precession. (It is this process 
that is responsible for the Earth's (forced) 'precession of the 
equinoxes' on a 26,000 year timescale, as the Sun and Moon 
exert a torque on the Earth's centrifugal deformation ). T he 



magnitude of the torque acting on the star is given by ( Sold- 
3tein 198C| ): 



T=^A/nsin2/?, 



(94) 



where jS is the angle that the angular momentum of this star 
makes with the normal to the plane in which the stars move, 
M is the mass of the other star and r denotes the stars' 
separation. This torque acts perpendicular to the plane con- 
taining the stars and J . We can write the quantity A In in 
terms of the rotation frequency Q. using equation ( |l3| ) ; 

This gives 



4^3 



(96) 
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As the two stars approach one another this torque will 
grow and change in orientation. If the stars pass very close to 
one another the steep r^'^ factor will give rise to an almost 
impulsive torque, acting when the stars are at and close 
to periastron, where the factor is at a maximum. If the 
interval in which the bulk of the angular momentum transfer 
takes place is much less than the spin period of the star the 
transfer will take place at nearly constant rotational phase, 
i.e. nearly constant reference plane orientation. It follows 
that the star would then be set into free precession with a 
wobble angle of order 



59 '. 



SJ I 



T 



Inserting T as given by equation then gives 

3 



7 :7 sm2P- — 

4 V d / /cruE 



(97) 



(98) 



where d denotes the separation at periastron. 

We therefore see that collisions capable of producing 
significant wobble angles {9 ~ 10~^) could conceivably oc- 
cur, but would require very close encounters, with a pe- 
riastron passage of no more than ten neutron star radii: 
d ~ lOR. Of course, during such close passages relativis- 
tic effects will be important, e.g. Lense-Thirring precession. 
We need not pursue these here. All that we require is an 
order-of-magnitude estimate of the collision cross section for 
a significant interaction to occur. Note that equation ( |98[ ) 
immediately rules out tidal torques during near collisions 
between neutron stars and main sequence stars as a pump- 
ing mechanism, due to the large {d J> 10® km) periastron 
separation. 

There will be an interaction between the neutron stars' 
magnetic moments also, as each dipole will tend to align with 
the field of the other. The torque on a dipole of moment 
m in a field B due to the other star is of order mB. The 
field B scales as r~'^. An argument analogous to the above 
gravitational one then leads to a wobble angle 



10-- (^) (^) 



100 Hz 



(99) 



where Bi and B2 denote the polar field strengths of the two 
stars. We have parameterised in terms of field strengths ap- 
propriate for a collision between a recycled and non-recycled 
neutron star. This is much smaller than the gravitationally- 
induced wobble angle and need not be considered further. 

Having established that near collisions could excite free 
precession we must now consider an event rate for such close 
passages. Of course, given that the event rate for encounters 
between a neutron star and a non-compact star was low, it is 
clear that the event rate for such close neutron star-neutron 
star encounters will be extremely low. In fact it is straight- 
forward to show that no such near-collisions will occur over 
a Hubble time, using a simple model. Suppose there are A'^ 
neutron stars in a globular cluster of size Rgc- Let Voo de- 
note their average velocity when far apart. Then in a unit 
time this population sweeps out an effective volume of or- 
der NAvoo, where A is a collision cross-section. Then the 
probability of a given neutron star colliding with another in 
this interval is of the order of this volume divided by the 
globular cluster volume, i.e. of order NAvoo/Rgc- As there 
are A'^ such stars the probability of any one of them collid- 



ing with any other is then N times this giving a collision 
rate Avoa/ Rgc- As there are approximately 200 globular 
clusters in the Galaxy we obtain a Galactic collision rate 



t^collisioii . 



J^collis 



200 



ngc 



(100) 



If gravitational attractions were neglected, the collision 
cross-section would be of order d^ ~ (100 km) However, 
gravitational focusing will increase the effective cross sec- 
tion as described in Verbunt & Hut (1987): 



A^d-' 



1 



2GMu 



.d 



(101) 



where Mtotai denotes the sum of the masses of the two stars. 
The second term on the right hand side describes the effects 
of gravitational focusing. For the case of interest it is the 
dominant factor. Parameterising we find an event rate 



11 -i/AV 
z^coiiision ~ 10 yr ( Yo^ / 



2.8 M, 



f d \ /lOkmN 
UOOkm/ V Wcx) / 




1 pc 

Roc 



(102) 



Such an event rate as this makes further comment unneces- 
sary, save to say that we will not consider stellar collisions 
any further. 



9 CONCLUSIONS 

This paper represents a systematic analysis of the detectabil- 
ity of gravitational waves from freely precessing neutron 
stars. It is based upon a model commonly employed to de- 
scribe the Earth's own motion. Explicitly, the neutron star 
has been modelled as an elastic shell with a fiuid core, whose 
angular amplitude of free precession (the wobble angle) is 
limited by its finite crustal lattice strength. It has been 
known for some time that neutron star structure may well al- 
low detectable gravitational wave signals at Earth, but this 
is the first study to attempt to identify particular astro- 
physical scenarios in which such precessional motion might 
be brought about and/or maintained. 

Broadly speaking, our findings were pessimistic. It 
proved impossible to find astrophysical pumping mecha- 
nisms capable of giving steady gravitational wave ampli- 
tudes detectable by an Advanced LIGO. This was because 
of the limiting effect of dissipation mechanisms internal to 
the star, even when dissipation strengths several orders-of- 
magnitude smaller than theoretically estimated values were 
assumed. 

Two qualifications are in order. Firstly, the above con- 
clusions were reached for stars with oblate deformations. 
In the physically less likely case of a star with a prolate 
deformation, the effect of internal dissipation would be to 
increase the wobble angle. Such a situation is certainly in- 
teresting from the gravitational wave point of view, although 
the dynamics of such a star, possibly involving crust crack- 
ing when the wobble angle exceeds a critical value, are far 
from clear. Secondly, most of the pumping mechanisms con- 
sidered in this paper involved an externally generated torque 
being exerted on the star. There exists another possibility. 
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where the symmetry axis of the deformation shifts due to a 
smooth plastic deformation of the crust. For instance, this 
deformation might be caused by accretion-induced temper- 
ature or composition asymmetries, in the manner described 
by Bildstcn (1998) for non-prcccssing triaxial stars. These 
two possible modifications to our model are currently under 
investigation, to see if there may yet prove to be a way of 
powering a long-lived, freely precessing gravitational wave 
source. 
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